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Numerous experiments showed that on cold compression graphite transforms into a new superhard 
and transparent allotrope. Several structures with difTerent topologies have been proposed for this 
phase. While experimental data are consistent with these models, the only way to solve this puzzle 
is to find which structure is kinetically easiest to form. Using state-of-the-art molecular-dynamics 
transition path sampling simulations, we investigate kinetic pathways of the pressure-induced trans- 
formation of graphite to various superhard candidate structures. Unlike hitherto applied methods 
for elucidating nature of superhard graphite, transition path sampling realistically models nucleation 
events necessary for physically meaningful transformation kinetics. We demonstrate that nucleation 
mechanism and kinetics lead to A/-carbon as the final product. W^-carbon, initially competitor to 
M-carbon, is ruled out by phase growth. Bct-C4 structure is not expected to be produced by cold 
compression due to less probable nucleation and higher barrier of formation. 



V_yomputational materials design is a central grand challenge of modern science. Recent progress in the develop- 
ment of theoretical methodologies for structure prediction, especially evolutionary algorithms [T] , has led to impor- 
tant discoveries confirmed by subsequent experiments In addition to thermodynamically stable solid phases, 
metastable modifications are worth exploring both computationally and experimentally because of their sometimes 
superior physical properties. However, among the infinite number of all possible metastable phases for each compound, 
only a handful are synthesizable for kinetic reasons. A metastable phase is synthesizable if it has the highest kinetic 
likelihood of formation at given conditions. 

The recently revived debate [H-^ on the nature of the metastable product of cold compression of graphite [TUI - ITS] 
represents a special case of interest. The transformation of graphite to the thermodynamically stable cubic diamond 
occurs at high pressures and temperatures (> 5 GPa and 1200-2800 K) and in the presence of catalysts. At ambient 
temperature, however, a new phase is formed above 15-19 GPa |10H13| and the transition is accompanied by clear 
changes in the physical properties of carbon [H [7] . 

Experiments have reported an increase in the electrical resistivity [T^l [T^HTSj and optical transmittance [TH [TS] . 
Raman spectra measurements up to 14 GPa indicate that the frequencies of the e'^^ rigid-layer shear mode at 44 
cm~^ and the in-plane mode at 1579 cm~^ increase sublinearly under pressure with initial pressure coefficients 
of 4.8 and 4.7 cm~^/GPa, respectively [14]. Changes indicating the formation of a new phase have been also observed 
in the x-ray diffraction pattern [131 ITBHTS] and in the near _ftr-edge spectroscopy [13] . More remarkable is the increase 
in hardness, as evidenced by the ability of the new phase to indent diamond anvils |13| . 

The difficulty to resolve the crystal structure of this mysterious superhard phase, which we call " superhard graphite" , 
from experiments has stimulated theoretical efforts [1|-17] . Hexagonal and cubic diamond have been ruled out because 
of the absence of their characteristic Raman bands at 1330 cm~^ [TH]. Amorphization was dismissed because of the 
persistence of Bragg peaks traceable to the original graphite pattern [U [T^ . 

Theoretically, candidate structures with distinct topologies featuring different patterns, combining odd (5 and 7) 
and even (4, 6, and 8) rings, have been proposed. A low-enthalpy monoclinic structure called M-carbon (space group 
C2/to, Z=16) was discovered using evolutionary algorithm USPEX [50] and was later identified with "superhard 
graphite" [J]. All carbon atoms in this structure are four-coordinate and form 5-membered and 7-membered rings 
when the structure is projected on (010). Later, a body-centered tetragonal structure called bct-C4 (space group 
//dmmm, Z=8) [5l[6l[2T] was proposed as another candidate. The four-coordinate carbon atoms form 4-membered 
and 8-membered rings on the (001) projection. Shortly afterwards, an orthorhombic candidate structure called W- 
carbon (space group Pnma, Z=16) was proposed [7]. This structure is very similar to that of M-carbon, featuring 
5-membered and 7-membered rings when the structure is projected on (010). Very recently, an orthorhombic structure, 
called 0CI6-II [Cmmm, Z=16), combining even-membered rings (4, 6, and 8) was predicted using particle-swarm 
optimization method [8] and metadynamics simulations [9]. Using the latter method, it was demonstrated that a 
large variety of superhard structures can be expected [e.g. 0CI6-I, mC12, and mC32) [3]. 

The simulated x-ray diffraction patterns of the predicted structures (M-carbon, FF-carbon, bct-C4, and 0CI6-II) 
are in satisfactory agreement with experimental data [ll[Sl[71|^. Their respective bulk moduli (and hardness [S^) are 
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392.6 (82.7), 391.8 (83.1), 393.4 (82.0), and 408.4 (84.4) GPa [S^. Band gaps calculated using the GMK method [53] 
are 5.00 |24J, 4.39 |7j, and 3.78 eV [5], for M-carbon, M^-carbon, and bct-C4, respectively. Ah candidate structures 
show no imaginary phonon frequencies and can explain the experimental observations of hardness and transparency 
of cold compressed graphite. 

Physical properties and energetics of these phases have been theoretically well studied [IHH US HB] , but the existing 
evidence is insufficient to decide which one was obtained in experiments [TDHT5] . In order to clearly identify the 
product of cold compression of graphite, we need to explore graphite transformation kinetics and compute energy 
barriers for transition regimes characterized by phase coexistence. 

Investigations of transition routes to superhard graphite [51 [7] using nudged elastic band method provided inconsis- 
tent results, especially inferring that graphite to cubic diamond has lower transition barrier than transformations to 
other three possible structures [7] . In addition, Ref . [7| reports a different sequence of energy barriers (cubic diamond 
< M-carbon < bct-C4) than Ref. [B] (bct-C4 < cubic diamond < M-carbon), despite the use of the same computational 
method. These conflicting results highlight the intrinsic complexity of a reliable elucidation of transition mechanisms 
and associated barriers. 

In a complex system, like a three dimensional solid, the lowest energy state of a system is not always the most relevant 
[27j . The viability of any candidate structure as a metastable product is strictly connected to the corresponding 
transformation pathway and the associated transition barriers or dynamical bottlenecks [TJ [551 US] • Transition barriers 
can be computed and bottlenecks identified with (stationary) saddle points representing transition states on the 
potential energy surface |30j . In complex systems with very rugged energy landscapes, saddle points cease to be 
characteristic points of the free energy barrier [31] . The latter barrier encompasses a large set of configurations and 
the associated relevant degrees of freedom are very difficult to anticipate [H 1301 [31] ■ To overcome this problem, we 
can use transition path sampling, a method based on collecting ensembles of dynamical trajectories [H [55H5T] . 

RESULTS 

In this work, we elucidate the nature of the metastable product of cold compression of graphite and compute the 
associated energy barrier from true dynamical pathways. The mechanistic details of this pressure-induced phase 
transformation are determined at the atomic level with full account for relevant events of nucleation. To achieve a 
reliable and detailed mechanistic picture of the transition we employ a strategy that combines isothermal-isobaric 
{NpT) molecular dynamics with transition path sampling (TPS) [55H51] (see Methods). This methodology has proven 
very effective in the simulation of activated processes [H [301 [35] with phase coexistence and growth [351 (Ml [3^ • 

How TPS works Transition path sampling is a generalization of Monte Carlo procedures in the space of dynamical 
trajectories connecting two states (graphite and metastable superhard phase in this case) separated by a high barrier 
in a rough energy landscape [32H34] . Therein, the relevance of a transition pathway is biased by path probability. The 
overall simulation approach is iterative and develops from an initial trajectory. By analogy with conventional MC 
simulations, the first step of TPS consists of equilibrating an initial pathway and gradually converging the trajectory 
regime to more probable regions. Accordingly, the first trajectory does not need to be a probable one or include 
mechanistic details of the real transition regime. 

In a situation characterized by the possibility of many theoretically viable candidate structures, TPS is a very helpful 
simulation strategy because it allows starting from a regime of very low probability or even a model constructed 
from a mapping between the limiting phases of the transition of interest. The convergence is ensured by Monte 
Carlo sampling of the space of trajectories. In the case of cold compression of graphite, we started from a pathway 
connecting graphite to cubic diamond generated by propagation of configurations obtained from a modeling approach 
based on transforming periodic nodal surfaces calculated from a short Fourier summation for each structure |37|, 138] 
(see Methods). 

Finding the real transition route The importance of nucleation events during the compression of graphite was 
pointed out for the transition to diamond at high pressures and temperatures [i39 . In the following, we investigate 
the kinetics of graphite transformation at ambient temperature due to intrinsic differences between nucleation and 
growth patterns of predicted superhard structures. 

The first set of TPS runs was started from trajectories connecting graphite to cubic diamond. Since diamond 
(cubic and hexagonal) was ruled out by experiments '19] and due to high activation barrier separating diamond from 
graphite, starting with this pathway will not only optimize the search in the path ensemble but will help converging 
to a more probable end point structure. The appearance of a new different phase instead of diamond is allowed by 
setting a fiexible order parameter that can accommodate four-coordinate structures (see Methods). 
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The set of transition pathways harvested in the course of TPS iterations shows a quick departure of the trajec- 
tory regime from the coUective motion inherited from the geometrical modehng. The graphite-diamond trajectories 
underwent three drastic shifts in the transition regime. 

First, after 15-20 iterations, the end point of the initial trajectory, graphite to cubic diamond, is replaced by a 
polytype intermediate between hexagonal and cubic diamond (Fig. [ij. The appearance of this polytype is due 
to opposite sliding directions of graphene layers giving rise to layers of hexagonal diamond. This evolution of the 
transition regime is similar to the results of a molecular dynamics investigation of the pressure-induced transformation 
path of graphite of diamond at high temperatures [40] . 

After 30-40 iterations, an inset of new structural patterns made of 5-membered and 7-membered rings appeared 
within 6-membered rings network (Fig. [2] (d)-(f)). These new patterns are still made of four-coordinate carbon 
atoms and interface very well with cubic diamond. They are characteristic of two structures suggested as metastable 
product for the cold compression of graphite, M-carbon and H^-carbon. The small size of the inset, however, does not 
yet allow to discriminate between these two models (Fig. [2] (f), dotted circle). 

Upon further sampling of the pathway ensemble, the inset of odd-membered rings grew larger and converted the 
remaining cubic diamond regions into a structure identical to that of M-carbon. The new transition regime connects 
graphite to M-carbon (Fig. |3|. Therein, the onset of the transformation is marked by the formation of interlayer single 
C— C bond (Fig. |3] (a)). The formation of the nucleus triggers a series of bond formation events along [001]graphite 
in a zigzag fashion (Fig. [sj (b)). The breaking of tt- interactions within graphitic layers by the emergence of the latter 
zigzag chain induces layer corrugation that facilitates the formation of 5-membered rings (Fig. [3) (c)-(d)). The latter 
rings pair up, and additional 7-membered rings are formed (Fig. [Sj (e)). 

The degrees of freedom induced by the mobility of graphene layers under pressure theoretically allow to consider 
different ways of graphene layers buckling. The survival of an initial regime or the appearance of a new one is connected 
to the way of nucleating the high pressure modification and the corresponding free energy barrier. TPS is biased to 
efficiently converge to the relevant set of the path ensemble which corresponds to the real transition pathway regime 
[29|, l32H35j and bypasses less-favored regions of the energy landscape corresponding to other mechanisms, graphite- 
bct-C4 and graphite- FF-carbon in this case. To illustrate this feature of TPS, we performed separate simulation runs 
to evaluate the energy barrier of graphite to bct-C4 and to 1/F-carbon transition routes and clarify the associated 
reconstruction mechanisms in comparison with the M-carbon route. 

Viability of graphite to W-carbon pathway Despite its structural and energetic closeness to the M-carbon 
phase [7], the arrangement of odd-membered rings corresponding to the orthorhombic candidate structure, M^-carbon, 
only appears in 20 % of trajectories collected during TPS simulations performed starting from graphite to diamond 
trajectories, over a pressure range from 15—20 GPa. 

In order to investigate differences in terms of transition mechanism and kinetics, we modeled a distinct set of TPS 
runs. The starting trajectory was taken from the previous runs that converged to graphite to M-carbon. Although 
nucleation events triggering the transition along the graphite to W-carbon path (Fig. |4j (a)-(b)) are very similar to 
that of graphite to M-carbon (Fig. [3j (a)-(d)), subsequent growth of W-carbon slightly differs from that of M-carbon. 
The formation of a C— C bond bridging two graphitic layers, in the case of W-carbon, induced the appearance of 
other bonds along [Q01]graphite as shown in Fig. p| Unlike in the transition mechanism to M-carbon, the zigzag 
chain propagates in a segmented fashion (Fig. |4|^b)j. The key step in the evolution of graphite toward M-carbon or 
W-carbon is the evolution of a single nucleation event, formation of an inter-layer C-C bond, either as infinite zigzag 
(Fig. [Sf^c)) or segmented chain (Fig. |4][b)). The breaking of the chain implies different progress of phase growth and 
results in a higher transition barrier (see below for transition kinetics). The average length of graphite to W-carbon 
transformation is 9 ps while graphite to M-carbon is 8 ps or less. This refiects slower transition rate for W-carbon. 

The effect of strain caused by the nucleus on the graphite lattice upon the direct conversion of graphite into diamond 
has been pointed out to be responsible of the formation of metastable hexagonal modification rather than cubic in 
agreement with experimental observations J39j. The nucleus formation along the graphite to W-carbon transition 
induced larger distortions of the graphite structure at the interface if compared to the path to M-carbon. This 
explains the more frequent appearance of the latter structure during the sampling. 

Viability of graphite to bct-C^ pathway As demonstrated above, the graphite to bct-C4 path is a less favored 
route. The set of sampling runs of such a trajectory ensemble quickly shift the transition regime toward the more 
favorable one, graphite- M-carbon. To prevent the disappearance of the bct-C4 structure, additional constraints are 
introduced on the order parameter to keep the transition regime as a pure graphite to bct-C4 transformation (see 
Methods). 

A representative trajectory of the latter path is depicted in Fig. [5j Unlike the graphite- M-carbon transition 
mechanism, the onset of the reconstruction is marked by the formation of two C— C bonds defining a C4 square unit 
(Fig. [5] (a)-(b)). The perturbation induced by the C4 nucleus helps propagating this motif and the bct-C4 phase 
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grows until complete reconstruction of graphite into the body-centered tetragonal structure. 

Kinetics of graphite cold compression In order to shed light on the kinetics of the transformation of graphite 
under high-pressure and ambient temperature, enthalpy variation of different simulated transition regimes starting 
from graphite have been computed and depicted in Fig. |6] The graphite to bct-C4 transition barrier (Fig. [6j star 
line) scores as the highest with 221 meV/atom. The graphite to cubic diamond transformation (Fig. |6j circle line) 
has a barrier of 200 meV/atom which completely rules out the latter transition route for graphite cold compression. 
Despite the similarities between transformation routes to W-carbon (Fig. [6j diamond line) and M-carbon (Fig. |6j 
square line), the latter is favored with a barrier equal to 176 meV/atom, lower than 194 meV/atom for the former 
transition. 

The role of graphite layers sliding In the nucleation process, the buckling of graphitic layers is realized upon 
their corrugation followed by formation of new C-C bonds (Fig.[I]-(b), Fig.[2]-(a), Fig.[3]-(b), Fig.|4]-(a), and Fig.[5]-(b)). 
This is facilitated by layer sliding in graphite that brings carbon atoms into an appropriate stacking sequence for C-C 
bonds formation along [001]. Just prior to nucleus formation, the applied pressure causes graphite layers to slide in a 
fashion that determines the nature of the high pressure modification to be formed. In Fig.[6| the variation of enthalpy 
indicates a different evolution for graphite even before the formation of any nucleus for all investigated transformations 
(graphite to cubic diamond, to bct-C4, to PF-carbon, and to M-carbon). For t ~ ~ 600 fs, graphite to bct-C4 shows 
a higher enthalpy kink if compared to the two other candidates, W- and M-carbon. This slight difference in enthalpy 
(40 meV/atom for bct-C4 and 30 meV/atom for W- and M-carbon) reflects two different layer sliding mechanisms as 
shown in Fig. [7| The graphite to bct-C4 transition requires an eclipsed arrangement of next-neighboring layer along 
[001] and changes the layers stacking sequence from ...AB... into ...AA... (Fig. [7j (a)-(b)) in order to trigger layers 
buckling. 

The transformations from graphite to 5-1-7 structures on the contrary require a reduced amount of displacements 
to induce graphitic layer buckling (Fig. [7j (a)-(d)). Therefore, transformation routes leading to W- and M-carbon are 
more favored than the path to bct-C4 in terms of layer sliding mechanism in graphite. Additionally, pathways to W- 
and M-carbon are less demanding in terms in-layer distortions, necessary to initiate nucleation. 

The stacking sequence and atomic displacements in the case of pathways to W- and M-carbon imply minimal in- 
layer distortions in order to initiate nucleation. The path to M-carbon is favored over the one to W-carbon because 
of different growth processes making the former structure more probable to form. 

The effect of simulation box size has been pointed out in large-scale molecular dynamics simulations of graphite to 
diamond transition [39] . If the simulation cell is too small, the derived transition mechanism does not feature in-layer 
distortions because in-layer stress implies artificial sliding of graphitic layers. The dimensionality of our simulation 
box fairly enables the observation of relevant nucleation events and allows to capture distortions around growing 
nucleus without inducing artificial layer sliding as shown in Fig[3]-(c) and Fig|4]-(b). Our results underline the role of 
distortions around the nucleus and how they determine the kinetics and the overall direction of the transformation. 

The very recently proposed candidate 0CI6-II deserves special notes. This orthorhombic structure can be described 
as an alternation of diamond and bct-C4 slabs along [100] (Fig. 3 in Supplementary Information). The existence of 
4-membered rings suggests mechanistic features of buckling in a fashion similar to graphite-bct-C4 transformation. 
To reconstruct graphite into a structure with a topology of rings fused as 4 -I- 6 + 8 pattern, atomic displacements 
in terms of layers sliding are expected to be larger than in the mechanism leading to M- and FF-carbon (Fig. 3 in 
Supplementary Information). On the experimental side, Zhao et al. [8] demonstrated an accordance of the calculated 
x-ray diffraction pattern of 0CI6-II with the experimental XRD data collected for cold compression of carbon nanotube 
bundles. The inconsistency with experimental data on "superhard graphite" has been pointed out [8 . In our TPS 
simulations, the initial graphite to diamond pathway evolved into the graphite to M-carbon pathway, bypassing the 
0CI6-II structure, in spite of its structural closeness to diamond. 

DISCUSSION 

Unlike synthesis from the gas phase, the high-pressure synthetic route of superhard graphite depends on the par- 
ticular nucleation history which would favor one pattern over another at the stage of phase growth within graphite. 
The viability of transformation routes is not determined by the overall stability of a particular candidate structure as 
an only parameter, but rather by activation energy of formation. Local events of C-C bond breaking and formation 
control the kinetics of the transition. In this work, an atomistic and kinetic picture of cold compression of graphite 
is presented using a methodology based on molecular-dynamics simulations combined with transition path sampling. 
It illustrates how transition path sampling can be used to assess synthesizability of metastable phases by optimizing 
transformation pathways. Unlike static models, we provide an unbiased description of different dynamical routes and 
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evaluation of the associated transition barriers. The final product of this first-order transformation is unequivocally 
identified as a metastable monoclinic structure, M-carbon. Although a competitor, FK-carbon is kinetically less favored 
and differences in transformation rates can be traced back to different patterns of phase nucleation and subsequent 
growth. We demonstrate that large distortions caused by the nucleation mechanism of W^-carbon make it less likely 
to form than M-carbon. Another earlier proposed structures, bct-C4, is ruled out because it did not appear during 
trajectory sampling under cold-compression conditions. 

It is worth mentioning that transformation pathways to "superhard graphite" are investigated as homogeneous 
nucleation processes. The presence of structural defects in the form of interlayer bonds can facilitate nucleation. 
Consequently, heterogeneous nucleation implies lower activation barriers. However, differences between homogeneous 
and heterogeneous nucleation do not alter the main features of the transformations revealed in this work. Phase growth 
kinetics around lattice distortions and structural defects are expected to be similar to that of homogeneous nucleation 
model. Intrinsic differences between different routes to "superhard graphite" are determined by seed growth fashion. 
In order to investigate the effect of defects on the kinetics and atomistic mechanism of the transition to " superhard 
graphite" , transition path sampling molecular dynamics simulations of heterogeneous nucleation will be carried out in 
the future. All in all, the investigation conducted in this work offers new insights into the understanding of materials 
synthesizability under high pressure, in a situation of competing phases and in absence of better resolved experimental 
data. 

METHODS 
Molecular Dynamics Simulations 

Born-Oppenheimer molecular dynamics simulations were performed in the A^pT ensemble 30J at temperature T=300 
K and pressures ranging from 15 to 20 GPa. Constant pressure and temperature were ensured by the Martyna- 
Tobias-Klein algorithm [41) allowing for anisotropic shape changes of the simulation box. Interatomic forces are 
computed within the framework of density functional tight binding approach [131 US] as implemented in the CP2K 
code |44j . For the accuracy of the set of DFTB parameters i45^ used, see Supplementary Information. The time 
propagation of the system was performed using the velocity Verlet algorithm [30j . To ensure good time-reversibility 
an integration timestep of 0.2 fs was used. The simulation box contained 256 carbon atoms. The starting configurations 
in all simulations correspond to graphitc-2II structure. The presence of defects and surface effects on the transition 
mechanism and the nature of the product of the transition is not investigated because the main focus is set on 
homogenous nucleation in graphite compressed at ambient temperature. 

Transition Path Sampling 

The sampling starts from an initial trajectory connecting the limiting phases. A new trajectory is generated 
by selecting a configuration from the existing one and slightly modifying the atomic momenta. The modifications 
Sp are applied to randomly chosen pairs of atoms according to: p^^"" = tf'^ + 5p{-t ^ - - T^j)| 

and it"^™ = — 6p{~T^ j — ~T^i)/\{~T^ j — keeping both momentum and angular momentum conserved. The 

resulting atomic momenta are rescaled by a factor of \J E'^^f^/ (J^i in order to keep the total kinetic 

energy conserved. The propagation of the modified configuration in both directions of time {-t,+t) generates a new 
trajectory. Repeating this step provides a set of trajectories and the successful ones are collected and analyzed. 
We initialized TPS with trajectories representing four different transformation regimes: graphite— cubic diamond, 
graphite— bct-C4, graphite— M-carbon, and graphite— W^-carbon. Subsequent equilibration of the path ensemble at 
300 K and pressures ranging from 15 to 20 GPa in two simulation runs for each transition regime resulted into 8 
independent TPS simulation runs. The number of trial shootings for each run exceeded 500 and the acceptance ratio 
was between 40-60 %. This production part of TPS procedure resulted in an ensemble of more than 200 successful 
pathways, between 8-10 ps long, of which ^70 are independent uncorrelated trajectories. 
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Order Parameter 

The average coordination number (CN) of carbon atoms is used as an order parameter to distinguish graphite 
(CN=3) from high pressure modifications made of four-coordinate atoms. Ah structures discussed in the manuscript 
(M-carbon, W^-carbon, bct-C4, hexagonal and cubic diamond) have an average CN=4 within the first nearest neighbor 
coordination sphere. In order to narrow down the path ensemble sampling we used second and third nearest neighbors 
to distinguish between different four-coordinate structures Clearly, apart from a quick detection of graphite and 
different candidate structures (bct-C4, M-carbon, and FK-carbon), this order parameter is capable of accommodating 
not only the three possible metastable phases but allows for starting transition pathways sampling from unlikely 
trajectories, like graphite to cubic or hexagonal diamond. The order parameter does not impose any bias on the 
evolution of the trajectory sampling. 

Topological Models 

We model transformation paths by transforming periodic nodal surfaces (PNS) [J [371 13H]- The reciprocal space 
approach implements short Fourier summations to define a family of surfaces according to the formula: f{x,y,z) = 
J2h k I \\Shki\\cos{2Tr{hx + ky + Iz) — cthki)- Shki is a geometric structure factor and ahki is the corresponding phase. 
The surface corresponding to f{x,y,z) = is called PNS. PNS of the transforming phases {A and B) are used 
as starting points for the definition of a geometric model for the transition. The choice of a common cell with 
a constant number of atoms and the periodicity of the model ensure commensurability of the two phases. After 
transformation of the reflections of the new setting of the cell, the transition can be formulated as a migration from 
one structure to the other along a coordinate s providing weighted linear mixing of the two functions: fAsix, y, z) — 
swAfA{x,y,z) + (1 — s)wb fB{x,y, z), s £ [0,1] To derive the transition state, configurations obtained from the 
geometric model in the range s G [0, 1] are propagated in molecular dynamics simulations at 300 K. Random velocities 
are assigned to atoms and the system is left free to evolve toward either A ov B without imposing any bias. Reversing 
the sign of the time coordinate allows reaching the other phase with a finite probability. 
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FIG. 1: Snapshots from a dynamical trajectory collected from transition path sampling connecting (a) graphite to (f) a polytype 
intermediate between cubic and hexagonal diamond. The buckling of graphite layers in initiated by the formation of C— C bonds 
along [001]graphite (b)-(c). Domains of cubic diamond are formed with different orientation (d)-(e). The interface between the 
latter domains defines a region of hexagonal diamond (e)-(f) 



CKXKX><X>0 CKDD-GC><D=C CKX><)C>=C0<) 




CMy><j><y>o ^ocHy>^>o ^oo^>o3^ 

GO-CO^XX) COOO-CCXB GCKX>a>CO 

CM30-C^<X>O C^CKXh^^O ^<5><X><X>0 
(a) (b) (c) 




FIG. 2: Snapshots from a representative trajectory illustrating the evolution of the (a) graphite to (f) cubic diamond transition 
regime. The mobility of graphitic layers during the reconstruction creates (d)-(f) an inset of 5- and 7-membered rings (dotted 
circle) within a 6-membered rings network. This inset interfaces well with cubic diamond and represents the seed of the 
metastable phase resulting from the cold compression of graphite. 
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FIG. 3: Snapshots of a representative trajectory of the stable regime corresponding to the cold compression of graphite. A single 
event of (a) bond formation (dotted circle) between graphitic layers triggers (b) a series of bond formation along [001] graphite 
in a zigzag fashion (dotted rectangle). The latter contacts facilitate the formation of (c)-(d) 5-membered rings causing the 
corrugation of graphitic layers and inducing the formation of (d)-(e) 7-membered rings. 
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FIG. 4: Snapshots taken from a representative graphite to VF-carbon transformation pathway. The buckling of graphite is 
initiated by (a)-(b) formation C— C contacts along [001]graphiie in the form of finite size zigzag chains. Each zigzag chain 
facilitates the corrugation of graphitic layers inducing the formation of (c)-(e) 5- and 7-membered rings transforming graphite 
into (f) K^-carbon. 
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FIG. 6: Enthalpy variation of different simulated transformations of graphite under pressure (15 GPa) and ambient temperature. 
The graphite to M-carbon transformation route (square line) indicates a lower energy barrier than the transition to VK-carbon 

(diamond line). The possibility of graphite transformation into bct-C4 structure (star line) on cold compression is ruled out 
because of higher barrier than the graphite to cubic diamond transition (circle line) . 
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(c) (d) 

FIG. 7: Top view of graphitic layers sliding fashions parallel to (001) plane. The (a) graphite to (b) bct-C4 transition requires 
an eclipsed arrangement of next-neighboring layer along [001] and changes the layers stacking sequence from ...AB.. into ...AA... 
in order to trigger layers buckling (only 3 graphitic layers are shown for better clarity). The (a) graphite to (d) M-carbon 
transformation implies small atomic displacements to induce (c) the onset of graphitic layers buckling. 



